Undular bore theory for the Gardner equation 
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We develop modulation theory for undular bores (dispersive shock waves) in the framework of 
the Gardner, or extended Korteweg-de Vries, equation, which is a generic mathematical model for 
weakly nonlinear and weakly dispersive wave propagation, when effects of higher order nonlinearity 
become important. Using a reduced version of the finite-gap integration method we derive the 
Gardner- Whitham modulation system in a Riemann invariant form and show that it can be mapped 
onto the well-known modulation system for the Korteweg-de Vries equation. The transformation 
between the two counterpart modulation systems is, however, not invertible. As a result, the study 
of the resolution of an initial discontinuity for the Gardner equation reveals a rich phenomenology 
of solutions which, along with the KdV type simple undular bores, include nonlinear trigonometric 
bores, solibores, rarefaction waves and composite solutions representing various combinations of the 
above structures. We construct full parametric maps of such solutions for both signs of the cubic 
nonlinear term in the Gardner equation. Our classification is supported by numerical simulations. 
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PACS numbers: 47.35.Fg, 47.55. Hd, 92.10.Hm 

I. INTRODUCTION 

The Gardner equation 

u t + 6uu x — 6au 2 u x + u xxx = , 



G) 



is a fundamental mathematical model for the descrip- 
tion of weakly nonlinear dispersive waves in the situ- 
ations when the higher order nonlinearity effects, de- 
scribed by the cubic term — 0>au 2 u x , become important. 
It first arose as an auxiliary mathematical tool in the 
derivation of the infinite set of local conservation laws 
of the Korteweg-de Vries (KdV) equation [l| but has 
been shown later to describe nonlinear wave effects in 
a number of physical contexts including plasma physics 
@j 0L stratified fluid flows [3] and quantum fluid dynam- 
ics [5j. One of the most important and best known ap- 
plications of the Gardner equation is the description of 
large-amplitude internal waves (see 0, 0, and refer- 
ences therein) . The coefficient a in (JXJ) can be positive or 
negative depending on the physical problem under con- 
sideration. In the context of internal waves, this depends 
on the stratification, see In the particular case of a 
two-layer fluid it is always positive [8(. 

When a — the Gardner equation ([1]) reduces to the 
KdV equation 



Using the change of variables 
1 

w = u , : 

2a 



(3) 



one transforms |1} to the modified KdV (mKdV) equa- 
tion 



Wt — 6aw w x i + w x 







(4) 



(note the change of boundary conditions at infinity) . 

The Gardner equation ([T]) is invariant with respect to 
the transformation: 



1 

a 



(5) 



u t + 6uu x + u x 



0. 



(2) 



which makes the existence of solutions of different po- 
larity (e.g. "bright" and "dark" solitons) possible for the 
same system, depending on the initial conditions. This is 
markedly different from the properties of the KdV equa- 
tion, which admits, for a given set of coefficients, solitary 
wave solutions of a fixed polarity, independently on the 
initial conditions. 

The soliton solutions of the Gardner equation for both 
signs of a are well-known and, along with the usual KdV- 
type "bright" and "dark" solitons, include table-top soli- 
tons, breathers, algebraic solitons and kinks (solibores). 
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Much less is known about dynamics of undular bores de- 
scribed by the Gardner equation. This problem is of 
significant theoretical and applied interest and is of par- 
ticular importance in oceanography, where undular bores 
play a key role in the evolution of the internal tide (see 

ii). 

Undular bores are nonlinear expanding wave trains 
connecting two different basic flow states and exhibit- 
ing solitary waves near one of the edges. They are usu- 
ally formed as a result of dispersive resolution of a shock 
or an initial discontinuity in fluid depth and/or velocity 
(see e.g. [l(J, (HI) or due to a resonant interaction of a 
fluid flow with localized topography (see e.g. [12], [l3|). 
Formation of undular bores (also often called dispersive 
shock waves) is a generic physical phenomenon which has 
been observed not only in classical fluids but also in colli- 
sionless plasmas, Bose-Einstein condensates and nonlin- 
ear optical media (see [3] and references therein). 

The analytical description of undular bores is usually 
made in the framework of the Whitham modulation the- 
ory (TBI . [iH ] in which the asymptotic solution for the bore 
is sought in the form of a slowly varying periodic so- 
lution of the governing dispersive equation. The slow 
evolution of the modulation parameters (such as mean 
value, amplitude, wavenumber etc.) is then governed 
by a hydrodynamic-type system of averaged equations, 
called the Whitham equations. The modulation descrip- 
tion of the KdV undular bore was first constructed in 
the celebrated paper by Gurevich and Pitaevskii [l?} and 
was later generalised to other dispersive systems both 
integrable (see e.g. fl8l. [Tpl) and references therein) and 
non-intcgrablc ( [n|, M\3, ) . 

The modulation system for the KdV equation can 
be represented in the Riemann invariant form [TEl [l6j 
which plays the key role in the Gurevich-Pitaevskii an- 
alytical construction of the KdV undular bore. For the 
mKdV equation (j4} (both defocusing (a > 0) and fo- 
cusing (a < 0) cases) the modulation system in Rie- 
mann invariants was derived in (23j using direct averag- 
ing of conservation laws and non-trivial algebraic manip- 
ulations leading to the diagonal structure. The spectral 
(finite-gap) approach to the derivation of the defocus- 
ing mKdV modulation system in the Riemann form was 
used in (24| , where the modulation solution was obtained 
for the undular bore resolving the "cubic" wave break- 
ing singularity. The undular bore theory for the focusing 
(a < 0) mKdV equation (j4]) was constructed in (25[. It 
was shown in [25| that, along with the KdV-type cnoidal 
undular bores, in which the elliptic modulus m varies to- 
gether with the wave amplitude a from m — 0, a = at 
the trailing edge to m = 1, a = a + at the leading edge, 
a + > being the amplitude of the lead solitary wave, the 
focusing mKdV equation supports another type of mod- 
ulated solutions, termed trigonometric bores, in which 
m = throughout the whole wave train but the ampli- 
tude a 7^ and vanishes only at the trailing edge. It was 
also shown in p5| that the trigonometric bore is usually 
realised as part of a composite solution: either a com- 



bined cnoidal-trigonomctric bore or a combination of a 
trigonometric bore and a simple rarefaction wave. Sim- 
ilar composite solutions were constructed in (2(| for the 
complex mKdV equation (which is related to the defocus- 
ing nonlinear Schrodinger (NLS), rather than the KdV, 
equation). We stress that neither trigonometric bores 
nor composite modulation solutions exist in the KdV and 
NLS modulation theories. These new patterns owe their 
existence to the fact that the mKdV modulation systems, 
unlike the KdV and NLS modulation systems, are nei- 
ther strictly hyperbolic nor genuinely nonlinear [2q |. Of 
course, the latter is not surprising if one remembers that 
the KdV and mKdV modulation systems are related by 
a non-invertible quadratic mapping (23j . a modulation 
counterpart of the Miura transformation. 

The Gardner equation, similar to the KdV and mKdV 
equations, is a completely integrable system, which im- 
plies that the Riemann invariants are in principle avail- 
able for the associated modulation system. However, we 
are not aware of any publications containing a consistent 
and complete derivation of the Gardner modulation sys- 
tem in Riemann invariant form (we note that some par- 
ticular results for the Riemann invariants via the map- 
ping between the KdV and Gardner spectral problems 
can be found in H?})- Consequently, the full theory of the 
Gardner undular bores has not been constructed. Some 
analytical progress has only been made for the case when 
the coefficient a is sufficiently small so that the Gardner 
equation can be asymptotically reduced to the KdV equa- 
tion via a near-identity transformation [28j . The undu- 
lar bore solutions in this case are qualitatively similar to 
their KdV equation counterparts provided initial discon- 
tinuity is not very large. An interesting phenomenology 
of the Gardner undular bore solutions, beyond the KdV 
par adig m, was revealed in the numerical simulations in 
j29L l30j . where the problem of the transcritical flow of 
a stratified fluid was considered in the framework of the 
forced defocusing Gardner equation for a broad range of 
values for a < and for the external forcing amplitude. 

In this paper we derive the modulation system for the 
Gardner equation in the Riemann invariant form and 
construct a full classification of the asymptotic (t ^> 1) 
solutions to the Gardner equation ([1]) with the initial 
conditions in the form of a step 

«m>={£; ^S: < 6 > 

We consider both signs of the coefficient a for the cubic 
nonlinear term. 

In the KdV equation theory, the resolution of the 
step ([6]) occurs via the generation of an undular bore if 
u~ > u + or a rarefaction wave if u~ < u + . For the 
Gardner equation we show that, due to the form of the 
nonlinear term in ([1]), the structure of the solutions to 
the initial value problem (QJ , (|6|) also depends on the po- 
sitions of the initial step parameters u + , u~ relative to 
the turning point u = 1 /2a of the characteristic velocity 
6u(l — au) of the dispcrsionlcss limit of the Gardner equa- 
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tion. The full classification encompasses 16 possible cases 
(8 for each sign of a). The wave patterns encountered 
include normal ( "bright" ) and reversed ( "dark" ) cnoidal 
undular bores, rarefaction waves, solibores (kinks), non- 
linear trigonometric bores and various combinations of 
the above patterns. 

The structure of the paper is as follows. In Section 
II we undertake the derivation of two families of peri- 
odic solutions of the Gardner equation ([l} correspond- 
ing to two signs of a. The solutions are derived in the 
"natural" parametrization by considering the travelling 
wave ansatz u = u(x — Vt) in (TTJ) and reducing it to an 
ordinary differential equation m| = Q{u), Q{u) being a 
polynomial of the fourth degree having (generally) four 
distinct roots ui < u 2 < U3 < u&, only three of which 
arc independent. The ordinary differential equation is 
then integrated in terms of Jacobi elliptic functions and 
the harmonic (m — > 0) and soliton (m — > 1) limits for 
both families solutions are then invesigated, m being the 
modulus of the elliptic solution. Section III is devoted 
to the derivation of the Whitham modulation equations 
in Riemann invariant form. For that, we take advan- 
tage of the reduced version of the finite-gap integration 
method [l~8| to derive the "spectral" representation of the 
periodic solutions obtained in the previous section. The 
outcome is the set of relationships between the spectral 
parameters ri,r 2 ,r3 and the dependent set u\,u 2 , 113,114 
characterising the periodic solution (two possible sets of 
relationships u(r) are derived for each sign of a — this is 
a consequence of the invariance of the Gardner equation 
with respect to the transformation (JSJ). The Whitham 
modulation equations are then derived for which rj 's are 
the Riemann invariants. In Section IV, based on the re- 
sults obtained in Sections II and III, we construct the full 
classifications of the solutions to the evolution of an ini- 
tial discontinuity problem for the Gardner equation with 
a > and a < 0. In Section V we draw conclusions 
from our analysis and outline possible applications of the 
obtained solutions. 



II. PERIODIC SOLUTION OF THE GARDNER 
EQUATION 



side of ([8]) as a "potential curve" for the nonlinear oscil- 
lator described by ©. Let Q(u) have four real roots 



Ul < U 2 < U 3 < W4 



(9) 



(the case when Q(u) has two real and two complex con- 
jugate roots corresponds to modulationally unstable so- 
lutions [23[ so we do not consider it here). The roots (O 
are obviously related by the condition 



2 

Q 



(10) 



and hence only three of them are independent. It is still 
convenient to keep all four Uj 's in the subsequent formu- 
lae to preserve the symmetry of the expressions. 

We should distinguish between two qualitatively dif- 
ferent cases. 

(a) Let a > 0. Then the periodic solution corresponds 
to the oscillations in the interval 



u 2 < u < u 3 , 
where the polynomial Q(u) is positive and 



(11) 



>/5(£-£o) = 



du 



lu yjiu — ui)(u - u 2 )(u 3 - u)(u4 - u) 

(12) 

The possible configurations of the "potential" curve Q(u) 
corresponding to qualitatively different travelling wave 
solutions are shown in Fig [1] 

The integral in (|12[) can be expressed in terms of the 
incomplete elliptic integral of the first kind and its inver- 
sion yields, after some algebra, the solution in terms of 
Jacobi elliptic functions: 



U 2 



(u 3 - M 2 )cn 2 (fl,TOi) 
^sn^rm) 



1 



U3 —U2 
U4 



where 



y/a(u 3 - ui){u A - u 2 )(x - Vt)/2 . 



(13) 



(14) 



We start with a direct derivation of the periodic trav- 
elling wave solution of the Gardner equation ([1]). Gen- 
eral expressions for such solutions can be found in (3l| . 
Here we need a more detailed description suitable for our 
subsequent development of the undular bore theory. In- 
troducing the substitution 



u = u(£), £ = x — Vt, 



(7) 



and integrating twice we arrive at a nonlinear oscillator 
equation 

ul = au i -2u 3 -Vu 2 + Au + B = Q{u), (8) 

where A and B are the integration constants. We shall 
sometimes refer to the polynomial Q(u) in the right-hand 



mi = 



(«3 - U 2 )(U 4 - Ml) 



(15) 



{Ui - u 2 )(u 3 - ui) 1 
and V is given by 

V = a(uiu 2 + U1U3 + M1U4 + u 2 u 3 + u 2 u 4 + U3W4). (16) 
The wavelength is given by the formula 

4K(mi) 



L 



^a(u 3 - iti)(u 4 - u 2 ) ' 



(17) 



where K(mi) is the complete elliptic integral of the first 
kind. The soliton limit toj — > 1 can be achieved in one 
of the two ways: when u± — > u 2 or when U3 — > U4. 
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FIG. 1: (Color online) Potential curve Q{u) configurations 
for the travelling wave solutions of the Gardner equation with 
a > 0. The oscillations occur between ui and U3. (a) peri- 
odic (elliptic) solution; (b) bright soliton; (c) dark soliton; (d) 
solibore; (e) linear wave. 



When U2 — > ui we obtain the "bright" soliton of eleva- 
tion propagating against a constant background u = U2 
(see Fig. lb), 



u 3 - Ul 



cosh 2 6 



sinli 

u 4 —Ui 



(18) 



If, further, one has Ui — 7/3 <C u 3 — «2, then the soliton 
(TT5|) becomes a wide, "table-top" soliton. 



Analogously, for u 3 
at £ = and obtain 

u(0 = u 4 



Ui we choose £0 so that it = u% 



Ui - u 2 



cosh 2 e 



■ sinh 2 6» 



(19) 



This is a "dark" soliton on the constant background u = 
Ui (see Fig. lc). If u 2 — Hi <C u 3 — «2 it assumes the form 
of a depression counterpart of the "table-top" soliton. 

If both U2 — > iti and U3 — > U4 then the polynomial 
Q(u) in the right-hand side of © has two double roots 
(see Fig. Id) , which implies that the solution assumes the 
form of a kink (a "solibore"). To study this limit, it is 
convenient to choose £ in such a way that u = (ui+Ui)/2 
at £ = 0. As a result, an elementary integration of © 
yields two possible solutions 



u(0 = Ui 



Ui — Ul 



exp[± v / a(u 4 — + 1 



(20) 



The lower sign corresponds to the kink with u — > it 4 as 
£ — > — 00 and u — > U\ at £ — > 00; the upper sign yields 
the "anti-kink" with u — » ui as £ — > — 00 and u — > Ui 
at £ — s- 00. As follows from (fT0|) . the limiting constant 
states ui and Ui are related by the condition ui + U4 — 
1/a. The speed of the kink (solibore) propagation in both 
cases is c = oT 1 + 2auiUi, which agrees with the shock 
speed obtained from the first conservation law ut + (3u 2 — 
2au 3 ) x = of the dispersionless limit of the Gardner 
equation. 

When U3 1*2 (mi — > 0) (sec Fig. lc) the cnoidal 
wave (|13p asymptotically transforms into a linear har- 
monic wave 

u = u-2 + \ {uz — U2) cos[k(x — Vt)], 

k = \J a(u 2 - ui)(ui - u 2 ), (21) 

V = 4u2 + a(uiUi — 3u 2 )- 

(b) Let now a < 0. Then periodic solution corresponds 
to the oscillations in one of the two intervals, 



ui < u < u 2 or u 3 < u < Ui, 



(22) 



where the polynomial Q(u) assumes positive values. The 
possible configurations of the potential curve Q(u) are 
shown in Fig. [5] 

First we consider the case 



Ul < u < U2, 



(23) 



so that 

VWks - Co 



du 



1 y/{u- ui)(u 2 - u)(u 3 - u)(u 4 - u) 

(24) 

A standard calculation yields 



U = 1*2 



(u 2 - Ui)cn 2 (6>, m 2 ) 



L sn z (c/, m,2) 



(25) 
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where 



VM(u 3 -ui)(u 4 -u 2 )(:r - Vt)/2, (26) 



m 2 



(Uj - u 3 )(u 2 - Ui) 
(u 4 - U 2 )(u 3 - Ml)' 

Now the wavelength is given by 

4K(m 2 ) 



y/\a\(u3 - Ui)(u4 - U 2 ) 

In the soliton limit u 3 — ¥ u 2 (m 2 — > 1) we get 

u 2 - u\ 



(27) 



(28) 



u = u 2 



cos\i 6 
V = 2u 2 + a(u\U4 — 3^2) 



2is^smh 2 (9' 

U4—U2 

2^ 



(29) 



This is a "dark" , depression soliton. 

The limit m 2 — > can be reached by two ways. 
(1) If U2 u± (see Fig. 2d) we get asymptotically 



u = u 2 — ^(u 2 — Mi) cos[fc(a; — Vt)], 

k = y|o|(u3 - "i)(w4 - ui), 
V = iux + a(u3ii4 — 3i* 2 ). 



(30) 



This is a small-amplitude harmonic limit. 

(2) If Ua = U3, but ui 7^ u 2 (see Fig. 2e) then we arrive 
at the nonlinear trigonometric solution 



u = u 2 



(u 2 — ui) cos 2 

1 + 2ia^"i sin 2 

M3-U2 



where 



= V|a|(u3-«i)(ua-«2)(a; - Vt)/2, 
V = 4u 3 + a(wiW2 — 3w 2 ). 



(31) 



(32) 



If we take the limit u 2 — ui -C i* 3 — ui in this solution, then 
we return to the particular case of the small amplitude 
limit pop with 114 = 113, so that k = \/\a\ (1*3 — u\) and 
V = 4u 3 + a(u\ - 3w§) = Am + a{u\ - 3u 2 ). On the 
other hand, if we take here the limit u 2 — > u 3 = 1*4, then 
the argument of trigonometric functions becomes small 
and we can approximate them by the first terms of their 
series expansions to get the dark algebraic soliton 



u = u 2 



u 2 — ui 



l + H(u 2 -«i) 2 (v-V^) 2 /4' (33) 
V = 2u 2 (l + am). 



Now we consider the case 

U3 < U < 1*4, 

so that 

>/H(£-£o) = 



(34) 



du 



^{u- ui)(u - u 2 ){u - u 3 )(u 4 - u) 

(35) 



Again, the standard calculation yields 

(114 - w 3 )cn 2 (#, 7712) 



u = u 3 + 



In the soliton limit u 3 — > u 2 (m — > 1) we get 

U4 — u 2 



where 



u = u 2 H q 

cosh 9 + sinh 2 6 

IL2—U1 



vVl( M 2 - "i)(u4 - u 2 )(x - Vt)/2, 



V = 4u 2 + a(u\U4 — 3u 2 ). 



(36) 



(37) 



(38) 



This is a "bright" elevation soliton. 

Again, there are two ways for getting the limit m 2 — > 0. 

(1) If U4 — > U3 (see Fig. 2f), then we obtain a small- 
amplitude harmonic wave 



where 



u = U3 + i(u 4 — U3) cos[fc(x — Vt)], 

k = y/\a\{u 3 - zti)(u 3 - u 2 )/2, 
V = AU3 + a{u\u 2 — 3it 3 ). 



(39) 



(40) 



(2) If u 2 — > u\ (see Fig. 2g), then we get another non- 
linear trigonometric solution 



u = u 3 + 



(114 — U3) cos 2 ( 
1 + i^a-sin 2 ! 



where 



6 = y/\a\(u 3 - ui)(u 4 -ui)(a; - Vt)/2, 
V = 4ui + a(w 3 it4 — 3it 2 ). 



(41) 



(42) 



If we assume here 114 — u 3 -C U4 — u%, then we reproduce 
the small amplitude asymptotics (|39j) with u 2 = ux, so 
that k = y/\a\(u3 — ui), V = 4m 3 + a(u\ — 3u 3 ). On the 
other hand, Eq. (|41l) in the limit u 3 — > u 2 = ui reduces 
to the algebraic bright soliton solution 



u = 1*1 + 



M4 — Ui 



l + H( U4 -wi) 2 (v-^) 2 /4' (43) 
V = 2iti(l + au 4 ). 

This completes the classification of stable periodic so- 
lutions and their limiting cases of the Gardner equation. 



III. SPECTRAL PARAMETRIZATION OF THE 
PERIODIC SOLUTION 

A. Motivation 

The periodic solution derived in the previous section is 
parameterized by four "integrals of motion" u\ < u 2 < 
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FIG. 2: (Color online) Potential curve Q{u) configurations for the travelling wave solutions of the Gardner equation with 
a < 0. (a) periodic (elliptic) solutions: ui < u < U2 or uz < u < 714; (b) bright soliton; (c) dark soliton; (d) linear harmonic 
wave, 7712 = 0, a = propagating about the background u = ui =112; (e) nonlinear trigonometric wave, 777.2 = 0, a 7^ 0, 
Mi < u < 712; (f) linear harmonic wave, 7712 = 0, a = propagating about u = 713 = 774; (g) nonlinear trigonometric wave, 
7772 = 0, a 7^ 0, 773 < u < 174. 



U3 < Ui, which are related by the condition (|10j) . In 
a strictly periodic solution these parameters Uj are con- 
stants, but in a modulated wave, which we are interested 
in, they become slow functions of space coordinate x and 
time t. Their evolution is then governed by the Whitham 
modulation equations (see (l5l.Tl6l fl8j) which can be ob- 
tained by averaging the conservation laws of the Gard- 
ner equation over the periodic solution family ([5]), and 
which, generally speaking, have the form of a quasilinear 
(hydrodynamic type) system 



dt 



2^ Vi i 



duj_ 
dx 



0, i, j = any three of {1,2,3,4}. 

(44) 

Here the matrix elements Vij are functions of u — 
(771,772,773,774) (note that one of the variables u% can be 
eliminated with the help of Eq. (fT0|) but then the sym- 
metry of the above expressions for the periodic solution 
will be lost). The modulation system in the form (|4"4")) 
would be, however, completely impractical due to the 
highly complicated structure of the matrix elements 
- this is already the case even for the KdV equation 
(see [HI, [ljl [HI ) . Fortunately, for the Gardner equation 



this system can be transformed to the Riemann diagonal 
form 



dr k dr k 



k = 1,2,3, 



(45) 



where r k , k = 1,2,3, are the Riemann invariants. This 
is possible due to the fact that the Gardner equation 
is a completely intcgrablc equation. Moreover, one can 
expect that, at least for a > 0, the Whitham system 
for the Gardner equation will be closely related (or even 
equivalent) to the Whitham system for the KdV equa- 
tion. Indeed, for a > the Gardner equation could be 
reduced, by a simple change of variables ([3]), to the defo- 
cusing mKdV equation, which, in its turn, is connected 
with the KdV equation by the Miura transform. As a 
result, for a > the mKdV- Whitham system in the Rie- 
mann form is equivalent to that of the KdV equation (the 
result first obtained in 1231) and the same is true for the 



Whitham- Gardner system [27[. The Whitham equations 
(|45|) can be readily solved analytically, providing the nec- 
essary modulation solutions. The problem, however, is 
that one still needs to know the dependence of Mi's on 
the Riemann invariants n, r2, r% for the Gardner equa- 
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tion to be able to find the modulations of the periodic 
traveling wave solutions obtained in the "natural" Uj- 
paramctrization. This dependence for a particular case 
of the travelling wave solution (TT5|) was found in [27J but 
the description in (2tJ , being merely an illustration of a 
more general theory, is too brief and somewhat incom- 
plete for our purposes, so below we present a detailed 
calculation, which also will not be restricted to the case 
a > 0. 

The transformation itj = Uj(r) is most conveniently 
found using the spectral theory of the Gardner equa- 
tion ([1]). The method of obtaining periodic (generally, 
quasiperiodic) solutions via the linear spectral problem 
associated with an integrable nonlinear dispersive equa- 
tion is usually referred to as the finite-gap integration 
method. It is based on the highly nontrivial properties 
of quasiperiodic solutions of soliton equations which have 
only a finite number of bands (gaps) in their spectrum 
when considered as potentials in the associated spectral 
problem (see e.g. j32|)- In the context of the Whitham 
modulation theory the advantage of the finite-gap inte- 
gration over the direct procedure of finding periodic so- 
lutions is that, if the endpoints of the spectral bands 
of the potential (quasiperiodic solution) are allowed to 
slowly vary with x and t, they become the Riemann in- 
variants of the modulation equations. The full finite-gap 
integration theory and the associated modulation theory, 
however, are quite technical and involve rather compli- 
cated algebraic-geometrical constructions on hypcrcllip- 
tic Riemann surfaces (see e.g. seminal paper |33j where 
this theory was developed for the first time for the KdV 
equation). However, in the single-phase periodic case of 
our interest a more simple, reduced version of the finite- 
gap integration method is available [l8j enabling one to 
derive the required Riemann invariant parametrization 
for periodic solutions of integrable equations associated 
with (2 x 2) linear spectral problems. 



B. Spectral theory 

The finite-gap integration method is based on the 
possibility to represent the Gardner equation (TTJ) 
compatibility condition of two linear systems (see, e.g. 

Ellll) 



where 
* = 



02 



u 



F G 
H -F 



V 



(46) 



A B 

C -A 



F = X, G=-(l-au), H 
-4A 3 - 2A(1 - au)u + u x , 



(47) 
(48) 



.4 

B = 4A 2 (1 - au) - 2Xau x - au xx + 2(1 - au) 2 u, (49) 
C = — 4A 2 u + 2Xu x — u xx — 2(1 — au)u 2 . 



This means that the condition ty xt = ty t x reduces to 
Eq. ©. 

Calculations become somewhat simpler if we transform 
this matrix form of equations to their scalar counterparts 
(see H|) 



4>xx = A-ip, ipt = -lB x tp + Bip x 



(50) 



where 



A 



B 



(A + g) 2 -(!-„»,„ -(£) (51) 



-4A 2 + 2A— - — -2(l-cmW (52) 
u u 



Then the second-order spectral equation in (|50)) has two 
basis solutions ip + and ip- and the "squared basis func- 
tion" 



9 = ^+^- 



(53) 



satisfies the third-order equation with a well-known inte- 
gral 



Ag 2 = P(X). 



(54) 



In the finite-gap integration method the periodic solu- 
tions are distinguished by the condition that P(X) is a 
polynomial in A. In our case we find that the one-phase 
periodic solution corresponds to the polynomial 



P(X) = - Y[(X 2 - A 2 ) = -(A 6 - Sl A 4 + S2 A 2 - s 3 ) (55) 
and 



.9 = A' - .91A + .92, 



(56) 



where the coefficients g\ and 172 are functions of x and 
t. Substitution of Eqs. ((55]) and ((56|) into Eq. ((54]) with 
A given by Eq. ((51"]) and equating of the coefficients of 
equal degrees of A at both sides of the resulting equa- 
tion yields a set of equations for gi, g2, as well as for u 
and its ^-derivatives. Elimination of g\, gi and of higher 
x-derivatives of u from these equations gives after a some- 
what tedious calculations the equation 

u\ = au 4 - 2u 3 + 2 Sl u 2 

^-^T v /(l-4aA 2 )(l-4aA 2 )(l-4aA 2 )] u 



1 

2^3 



1 - 2as 1 - 2a 2 (si - 4s 2 



Tx/(l-4aA 2 )(l-4aA 2 )(l-4aA 2 ) 



(57) 



The travelling periodic solution of Eq. (fT]) is then ob- 
tained by the replacement x n- x — Vt, where V = — 2s\ 1 
that is 

y = 2(A 2 + A 2 + A 2 ) = 2(r 1 +r 2 +r 3 ), (58) 
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where we have introduced r, = Af, i = 1,2,3. The pa- 
rameters Tj are expected to become the Riemann invari- 
ants of the Whitham modulation equations. 

We note that the representation ([57]) of the ordinary 
differential equation for a periodic solution, unlike that 
given by its equivalent ((8]), contains two possible signs. 
As a result, one needs to use different sets of relationships 
between {ui} and {r;} for different types of solutions, i.e. 
the mapping {ui} n- {r^} is not one-to-one. To express 
the original parameters Uj (the zeroes of the polynomial 
in the right-hand side of equation ©) in terms of the 
spectral parameters r; = Af we compare the two forms 
of the same ordinary differential equation defining the 
the periodic solution, namely, Eq. (JSj) and Eq. (f57]) . 

First let us consider the case a > 0. The example 
of the mKdV equation (see (24[) suggests the following 
expressions which can be verified by direct calculations: 

1 



Mi = — (l — \/l — 4oti — y/1 — Aar-y + y/1 — Aar 3 ) 
2a 

1 

277 
1 

2ct 
1 



u-2 = — (l — vT — 4oti + a/1 — 4ar 2 — V 1 — 4ar 3 ) 
2a 



U3 = — (l + VT— 4qt7 — yi- 4ar7 — V 1 — 4ar 3 ) 
2a 



u 4 = — (l + Vl - Aari + Vl - 4ar 2 + \/l - 4ar 3 ) 



(59) 



in the case of the upper sign in Eq. (|57p and 



It! = 


1 

277 


(1- 


"2 = 


1 

277 


(1- 


U3 = 


1 

277 


(1 + 


U4 


1 

2a" 


(1 + 



(60) 

in the case of the lower sign. In both cases the zeroes Ui 
are ordered according to © provided 

A? < A 2 < A3 or, equivalently, n < r 2 < r 3 . (61) 

For both cases (|5"9"|) and (|60[) the inverse formulae are 
simply 



n = — («i + « 2 )(u 3 + 04), 

^2 = + « 3 )(u 2 + u 4 ), 

r3 = ^(u 2 + u 3 ){ui + m). 



(62) 



The existence of two sets ([59^1 . (JHUJ) of the travelling wave 
parameters corresponding to the same set of the spectral 
parameters rj is due to the invariance of the Gardner 
equation with respect to the transformation ([5])- Indeed, 
the set of relationships (JBDJ) can be obtained from (f5"9")l 



by applying the transformation Uj — > 1/a — Uj and then 
reordering the resulting set. 

In the case of a < the expressions remain the same 
but their order corresponding to © is different. For the 
upper sign in (|5T[) we obtain 





1 


Ul = 




277 




1 


U-2 = 




277 




1 


U3 = 


2a 




1 


U4 = 


2a 



(63) 



and for the lower sign 





1 


Ul = 




277 




1 


U-2 = 




2a 




1 


U3 = 


277 




1 


U4 = 






277 



(64) 



Now the inverse formulae arc 



(65) 



r\ = —{u\ + u 4 )(u 2 + u 3 ), 
ri = —(ui + u 3 )(u 2 + Ui), 

It is essential that the expressions (fT5|) and (|27l) for the 
modulus of the elliptic function for the cases a > and 
a < respectively, reduce to the same formula in terms 
of r fe 's: 



mi = m-2 = m = 



r 2 - ri 



(66) 

and in both cases the wavelength is given by the formula 

L=4SL (67) 

The periodic solutions obtained in the previous sub- 
section can now be written down directly in terms of 
rj's. This would lead, however, to rather cumber- 
some expressions so it is better to keep the original Uj- 
paramctrization in the periodic solutions and use the re- 
lationships (|59p . or (j6"3")l , for imposing slow mod- 
ulation more conveniently represented in terms of Tj's. 



C. The Whitham modulation equations 

The Whitham modulation equations in the Riemann 
form (|45|) can be derived using the well-established pro- 
cedure of averagin g th e generating equation for conserva- 
tion laws (see e.g. [la|)- This procedure for the Gardner 
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equation, however, is not as straightforward as it is for 
the KdV or mKdV equations, so, to avoid lengthy calcu- 
lations, we make a plausible assumption that the roots 
rj of the "spectral polynomial" (f55j) are the Riemann in- 
variants of the associated modulation system (|45|) as it is 
the case for the related KdV and mKdV equations and 
other integrable systems. Then we observe that expres- 
sions (|59")) -(1 M| ) do agree with the corresponding KdV (l5j 
and mKdV ([23[) expressions in the limits as a — > and 
a — > ±00 respectively. We also note that the particular 
set (|59")l of the relations between Ui's and r,'s was actu- 
ally obtained (up to some obvious misprints) in [2tJ using 
the traditional finite-gap method; the remaining expres- 
sions (f6T)|) . (|M|) have the same structure and can be 
simply derived as extensions of that result. 

Having established the modulation Riemann invari- 
ants, the expressions for the corresponding characteristic 
speeds Vi(r) can be derived directly from the wavenum- 
ber conservation law, by-passing thus the detailed averag- 
ing procedure. To this end we consider the wavenumber 
conservation l aw, which is a generic modulation equation 
(see, e.g., EHEl), 



k t + (kV) a 



0. 



(68) 



where fc(r) = 2ir/L and V(r) are the wavenumber and 
the phase velocity expressed in terms of the Riemann 
invariants. Since Eq. (|68[) must be consistent with the 
diagonal system (|4"5)l , one readily obtains the "potential" 
representation 



1 



L 



o l v, a 



_d_ 

dn 



(69) 



The function V(r) is given by (|58[) and the dependence 
L(r) by (J6TJ) . As a result, Eqs. (j69j) yield the Whitham 
characteristic velocities 



2(ri+r 2 + r 3 ) + 



v 2 = 2(n +r 2 + r 3 ) - 
v 3 = 2(n + r 2 + r 3 ) + 



4(r 2 - ri)K(m) 

E(m) - K(m) ' 
4(r 2 — ^i)(l — m)K(m) 

E(m) - (1 - m)K(m) ' 
4(r 3 - r 2 )K(m) 
E(m) 



(70) 



One can see that the characteristic velocities ([70)1 coin- 
cide with those for the KdV modulation system (see e.g. 
[l~8l HH). This, however, does not imply that the dy- 
namics of the modulated periodic waves in the Gardner 
equation will necessarily be the same or even qualita- 
tively similar to the KdV case. Indeed, to obtain the 
modulated solution for the Gardner equation, one needs 
first to convert the solution for ri(x,t) into the dynam- 
ics of the original modulation parameters Ui(x, t) and 
then substitute u^'s into the relevant periodic solution 
from Section II (which could be drastically different from 
the KdV cnoidal wave). Moreover, since the mapping 
{ffc} l— ^ {ui} is not one-to-one, the same modulation so- 
lution {rfc(x, t)} can give rise to two completely different 



modulations {ui(x,t)} of the periodic solutions to the 
Gardner equation. This becomes evident already on the 
level of the dispersionless limit of the Gardner equation, 



u t + 6u(l — au)u x = . 



(71) 



which is related to the dispersionless limit of the KdV 
equation (the Hopf equation) 







(72) 



± of the quadratic equation 
= 1/a). These two con- 



via the quadratic mapping u ^ w specified by the func- 
tion w = u(l — au). Indeed, a constant solution w = a 
of equation (|T2"|) gives rise to two different constant solu 
tions of (fTTj) found as roots u 
w(u) = a (obviously u + + u 
stant solutions can also be combined into a weak solution 
in the form of a propagating step: u = it_ for x < 6at 
and u = u+ for x > 6at (we note that this step solution 
transforms into a smooth kink (or solibore) solution ([20)) 
of the full Gardner equation if dispersion is taken into 
account). 

It is clear that a one-to-one correspondence between 
the solutions of (fTTj) and ([72^1 is only possible in one of 
the restricted domains of the function w(u): u < l/2a 
or u > l/2a, where w{u) is monotone. Thus, one can 
expect significantly different, compared to the KdV case, 
dynamics if the range of the initial function u{x 1 0) would 
include an open interval containing the turning point u = 
l/2a of the characteristic velocity of the dispersionless 
Gardner equation (|7ip . 

The above non-uniqueness in the correspondence be- 
tween the modulation solutions of the KdV and Gardner 
equations is due to the invariance of the Gardner equa- 
tion with respect to the transformation ([5])- As a result, 
the Whitham-Gardncr system in "natural" modulation 
variables, unlike the KdV- Whitham system, is neither 
genuinely nonlinear nor strictly hyperbolic (see [26| for 



the detailed analysis of a similar issue in the context of 
the closely related complex mKdV system). This results 
in the occurrence of much richer modulation dynamics 
for the Gardner equation than those for the KdV equa- 
tion. A very similar issue was also recently discussed in 
[23 | where the dam-break and lock-exchange flows were 
studied for the Miyatta-Choi-Camassa (MMC) system 
[3l| describing fully nonlinear long dispersive interfacial 
waves in a two-layer fluid. This is not surprising, of 
course, as the Gardner equation with a > represents 
a weakly nonlinear approximation of the MMC system 
obtained under an additional requirement that the layers 
depth ratio is close to the critical value defined by the 
square root of the ratio of the respective fluid densities. 

In conclusion we note that the Whitham modulation 
system associated with the Gardner equation with a < 
can be elliptic (this is possible since the spectral eigenval- 
ues Xi can be pure imaginary so the the squared Riemann 
invariants are negative and the mapping {r;} t— > {ui} 
is generally complex). It is indeed the case for the related 



focusing mKdV equation u t + 6u u x 



(see e.g. 
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[23( |) so modulational instability is generally to be ex- 
pected in this case. However, it was shown in (36j that 
for the focusing mKdV equation with real initial data, the 
modulation equations arising in the zero-dispersion limit 
of the initial-value problem, are hyperbolic which guar- 
antees modulational stability for such problems. This 
property was recently used in (25j to construct undular 
bore solutions to the focusing mKdV equation. Since the 
focusing mKdV equation and the Gardner equation with 
a < are related by the simple transformation ([3]), we 
shall be assuming hyperbolicity of the Gardner modula- 
tions for the resolution of a step problem considered in 
the next section. 



IV. CLASSIFICATION OF THE SOLUTIONS 
FOR THE STEP PROBLEM 

We now consider the Gardner equation ([1]) with the 
initial conditions in the form of a step, 

«m) = {£; III ^ 

It is clear that, due to the form of the nonlinear term in 
the Gardner equation, the structure of the solutions to 
the initial value problem ((T|), ([73"]) will strongly depend 
on the positions of the initial step parameters u + , u~ rel- 
ative to the turning point u = l/2a of the dispersionless 
characteristic velocity 6u(l — au). 

A. Key ingredients 

We first describe several particular solutions of the 
Gardner equation playing the role of "building blocks" 
in the full solutions to the dispersive Riemann problem 
®, (|73|) for different values of u ± . These solutions are: 
cnoidal undular bores, rarefaction waves, solibores (for 
a > 0) and trigonometric undular bores (for a < 0). 

a) Cnoidal undular bores: Gurevich-Pitaevskii solu- 
tion. 

The local structure of the simple undular bore is de- 
scribed by one of the periodic solutions obtained in Sec- 
tion II: solution (fT3|) for a > and solution (|25|) for 
a < 0. The corresponding modulations are expressed in 
terms of the parameters r\ , r 2 , r% satisfying the Whitham 
equations (|4"5)l . (fTUj) . It is clear that, in the problem of 
dispersive resolution of an initial discontinuity we are in- 
terested in the similarity modulation solutions where the 
modulation variables depend on s = x/t alone (both ini- 
tial data and the modulation equations are invariant with 
respect to the scaling transformation x — > Cx, t — > Ct). 

The classical Gurevich-Pitaevskii similarity solution of 
the modulation system (|4"B")) . ((70)) has the form (see [I?], 

Ell) 

ri = r + , r 3 = r~ , (74) 



where r + ,r are some constants, while the dependence 
r 2 (x,t) is given implicitly by 

x 

v 2 (r + ,r 2 ,r~) = - . (75) 

It is essential that, since r^ > r\, one must have 

r~ >r + . (76) 
A typical modulation solution is presented in Fig. 6. 



.!5 




x 



FIG. 3: (Color online) Typical behaviour of the Riemann 
invariants rj in the modulation solution for a simple cnoidal 
undular bore. 

The undular bore described by the Gurevich-Pitaevskii 
modulation solution (|74[) , (|75[) occupies an expanding re- 
gion x~ < x < x + whose edges x^ = s t propagate with 
constant velocities s . The trailing (harmonic) edge is 
defined by the condition m = (i.e. r 2 = r\ = r + 
see (f6"T|) and the leading (soliton) edge — by the condition 
m = 1 (i.e. r-2 = ^3 = r~). Then the velocities are 
found from {73}, dZOj) as 

s~ = v 2 \ r2=ri = 12r x - 6r 3 = 12r+ - 6r" , (77) 

s+ = v 2 \ r2=ra = 2n + 4r 3 = 2r+ + 4r" . (78) 

The above solution (|74|) . ([75]) coincides with the modu- 
lations in the undular bore arising as a result of the res- 
olution of an initial discontinuity: r(x > 0, 0) = r + and 
r(x < 0, 0) = r~ for the KdV equation r t + 6rr x + r xxx = 
[13]. However, in the KdV context the modulation so- 
lution ([74"]) . (fT5"j) uniquely characterises the asymptotic 
solution for KdV undular bore due to the one-to-one cor- 
respondence between the Riemann invariants {rj} and 
the physical parameters {uj} of the travelling wave solu- 
tion [13| . In the case of the Gardner equation ([T]) there 
are two possibilities for each signs of a due to differ- 
ent possible relationships between the {r^} and {ui} for 
different types of travelling wave solutions described in 
Section IIIB — see relationships (|5T))) , (|M| for a > and 
(f6"3")> . (fM)) for a < 0. The actual choice depends on the 
positions of the initial step parameters u + and u~ rel- 
ative to the turning point u = l/(2a) of the function 
w(u) = u(l — au) (see Section III C). We shall consider 
all possible cases in the next two sections. 
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b) Rarefaction waves. 

The rarefaction waves are asymptotically described by 
the similarity solution 



.r 

7 



6u(l — cm) 



(79) 



of the dispcrsionlcss limit (|71l) of the Gardner equation. 
The two possible roots of (|79|) describe two types of rar- 
efaction waves: "normal" and "reverse". We first con- 
sider the case a > 0. 

The normal simple rarefaction wave connects two con- 
stant states u = u l (left) and u = u r (right) satisfying the 
condition ui < u r < I /2a and has a structure similar to 
the rarefaction solution of the "rightward-propagating" 
KdV equation, in which du/dx > 0. It is described by 
the equation 



u = u l for x < s l t . 



1 



1 



1 - 



2ax 



2a \ V M 

u r , for x > s r t . 



for s t < x < s r t . 



(80) 



We shall be using a symbolic diagram {u l RW — » u r } 
for the normal rarefaction wave. 

The reversed simple rarefaction wave connects two con- 
stant states u = u l and u — u r satisfying the condition 
u r < ui < I /2a and is described by the equation 



u = u for x < sit . 



1 



1 



1 



2ax 



2a \ V 3i 

it r , for x > s r t . 



for s t < x < s r t , 



(81) 



The reversed rarefaction wave is similar to the rarefaction 
wave in the "leftwards-propagating" KdV equation, in 
which du/dx < 0. The symbolic diagram for this wave 
is {ui 4r- RW u r }. 

The speeds s l,r of the left and right boundaries for both 
normal and reversed rarefaction wave are given by 



6u Lr (l 



au 



s l < s T 



(82) 



Solutions ([8Tj and (|80|) have weak discontinuities at the 
corners x = s l ' r t. These are smoothed out by small- 
amplitude oscillatory wavctrains which are generated if 
the dispersive term of the Gardner equation is taken into 
account (cf. HH for the KdV case). 

For a < 0, the rarefaction waves are described by the 
same solutions, however, formula ([50)1 would then de- 
scribe the reversed wave with u r < u l < — l/(2\a\) and 
(|5T|) — the normal one with u r > u l > — l/(2|a|). 

c) Solibores (a > 0). 

The solibore (kink) solutions for the Gardner equation 
with a > are given by formulae (f2"U|) . Solibores pro- 
vide the smooth transition between two constant states 
ui and u r satisfying the condition ui + u r = l/a. Using 



the terminology introduced above for rarefaction waves, 
we shall refer to the solibore as "normal" when we have 
du/dx < ("-" sign in solution and "reverse", when 
du/dx > ("+" sign in (|2T))) ). The corresponding dia- 
grams are: {ui SB — > u r } for the normal solibore and 
{ui SB Ur} for the the reversed one. 



d) Trigonometric undular bores (a < 0). 

This type of undular bores, not encountered in the 
KdV theory, was first reported in (25j where the evolu- 
tion of a step problem was studied for the focusing mKdV 
equation (see also a similar solution for the complex mod- 
ified mKdV equation in 26] ) . The trigonometric undular 
bores of the Gardner equation with a < are described 
by the modulated finite-amplitude nonlinear periodic so- 
lutions (pTtj) or (HH so that m = throughout the wave- 
train. At one of the edges of the trigonometric bore 
the amplitude vanishes and at the opposite edge it as- 
sumes some finite value. Generically, as will be explained 
later trigonometric undular bores are realised as parts of 
composite solutions (either a combination of cnoidal and 
trigonometric bores or a combination of a trigonometric 
bore and a rarefaction wave). As with other wave pat- 
terns arising for the Gardner equation, one can have two 



types of trigonometric bores: normal, {u TB 



and reversed, {u l TB u r }. The normal trigonometric 
bore is locally described by solution (pTTj) while for the 
reversed one solution (13TT) should be used. 
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FIG. 4: (Color online) The sketch of the "square" Rie- 
mann invariants n, ri, behaviour in the composite cnoidal- 
trigonometric undular bore 



A typical configuration of the Riemann invariants rj 
in the composite cnoidal-trigonometric bore is shown in 
Fig. 4. Since in the region of the trigonometric bore one 
has r2 = T\ = r, the corresponding similarity modula- 
tion solution v-zir, r, r~) = x/t is degenerate in the sense 
that it does not allow one to reconstruct uniquely the 
modulations for Uj(x, t) in the trigonometric bore and to 
provide the necessary matching between u l and u r across 
the bore. Therefore, the "square" Riemann invariants rj 
(f6"2"|) used so far, are not suitable for the description of 
trigonometric bores. Instead, motivated by the results in 
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[25j we introduce the classical Whitham combinations and so from (|83[) we get 



Ri = s(«2 + «a) = 2~ (1 + VI - 4arij 



#2 = |(«2 + u 4 ) = 2~ ( x - Vl-4ar 2 J 



(83) 



#3 = 5(^3 + u 4 ) = ^ (1 - VI - 4ar 3 ) 



(i?3 > R2 > Ri) for the normal trigonometric undular 
bore (see the relationships ([51])) and 



#3 = + «a) = ^- (1 + Vl + 4ar 3 ) , 
Za 



R2 = \{ui + u 3 ) = — (1 + VI + 4ar 2 ) , (84) 



#1 = 5 ("2 + «s) = 5^ (1 - VI -4an) , 

(i?i > i?2 > -R3) for the reversed trigonometric undular 
bore (see the relationships 

Obviously, the quantities Rj(rj) are the Ricmann in- 
variants of the modulation system (any function of the 
Riemann invariant alone is also a Riemann invariant). 

(i) Normal trigonometric bores. 

We now construct the modulation solution for the nor- 
mal trigonometric bore, where the oscillations occur be- 
tween the roots U3 and 1*4 of the travelling wave solution 
polynomial Q(u) ([8]). The amplitude in such a bore is 



a = Ua — U3 = 2(i?i - R 2 ) 



(85) 



and it gradually increases from a — at the left (har- 
monic) edge, say x = ar, to some nonzero value a = a r 
at the right edge x = x r . 

When a = we have from (|4"Tj) u = 1x3 = 1*4. Hence, 
since the trigonometric bore must match with u = u l at 
x = x l we obtain from ([531 

Ri = U4 = u l at ie = ar, 



(86) 



and, therefore, -R3 = u l everywhere within the trigono- 
metric bore. Hence, from (1831) we obtain 



(87) 



r 3 = u l (l - au l ) 



Now, within the trigonometric undular bore we have 
r% = r% = r(x/t), which is determined by the degenerate 
similarity solution of the Whitham system (|45[). 



(>2 



= 12r - 6r 3 = - 



Substituting the value of r 3 (|87[) we obtain 
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+ 6u l (l-au l ) 



(89) 



Ri 



Ri 




2a 2 



1 

2o" 



lax 
~3T 



2a 2 



1 

2o 



2ax\ 



(90) 



One can see that R\ + R% = 1/a for all a; in the trigono- 
metric bore. 

If a trigonometric undular bore is fully realised (i.e. 
is not part of the composite cnoidal-trigonometric bore) 
then at the leading edge x = x r it must assume the lim- 
iting waveform of a bright algebraic soliton (|43[) (other- 
wise the matching with constant or smooth external so- 
lution would not be possible) . This implies that we have 
M3 = U2 = ui and so R2 = R3 = u l at the leading edge 
x = x r . The algebraic soliton rides on the background 
ui so the relevant matching condition becomes 



Ul 



at x 



(91) 



which, by ui = U2 = U3 and the first formula in (j83[) 
implies Ri = u r at x — x r . Therefore, the trigonometric 
bore can only connect the states u l and u r satisfying the 
condition 



1 

* 

Q 



(92) 



i.e. a single isolated trigonometric bore can be realised 
as a result of the step evolution only in the special cases 
when the parameters of the initial step satisfy the condi- 
tion (j92[) (note that for normal trigonometric bore one 
must have u l > u r which follows from the ordering 
i?3 > Ri). Thus, for the Gardner equation with a < 
the trigonometric undular bores play the role similar to 
that played by solibores in the step problem for the the 
Gardner equation with a > 0. 

The speed s l of the trailing edge of the trigonometric 
bore is found from the condition that at the trailing edge 
R 2 = Ri = J- which implies by 



s = 6cra (1 — au ) . 

a 



(93) 



At the leading edge x = s r t of the trigonometric undular 
bore we have R2 = R3 = u ■ Then, substituting R2 = u , 
x/t = s r into the second equation (j90[) we obtain the 
speed of the leading edge 



6u l (l — au 1 ) = 6it' (1 — au r ) . 



(94) 



which coincides with the characteristic speed of the dis- 
persionless Gardner equation at u = u r . This implies, 
in particular, that the trigonometric undular bore can be 
joined at the leading edge to a simple rarefaction wave 
solution. 
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The amplitude of the algebraic soliton (|43|) at the lead- 
ing edge follows from ([8"5"j) where we set R 2 = R3: 

a r = 2(R 3 -R 1 )\ x = x r = 2{u l -u r ) = 2(2u l -l/a). (95) 

(ii) Reversed trigonometric bores. 

In the reversed trigonometric bore the oscillations oc- 
cur between the roots U\ and u 2 of the polynomial Q(u) 
in ©. The amplitude is given by 



u 2 -u 1 = 2{Ri - R 2 ) 



(96) 



In terms of the Ricmann invariants rj the modulation so- 
lution for the reversed bore is given by the same formula 
([89]). which is then translated to R/s d84j) as fcf. ([90l) ) 




1 / 2ax\ 

2 [ ~~sTJ 



(97) 



Similar to the normal trigonometric bore, the reversed 
trigonometric bore has a restriction (|92|) for the admis- 
sible boundary values u l and u b . At the leading edge 
x = x r the reversed bore assumes the limiting form of 
a dark algebraic soliton (|33|) . which has the amplitude 
a r = 2(2u l — I /a). The trailing and the leading edge 
speeds are given by the same expressions (|93]l and (|94|) 
respectively. 

As was already mentioned, the trigonometric undu- 
lar bore (normal or reversed) can occur as part of the 
composite cnoidal-trigonometric bore. In that case it is 
realised only partially and does not contain an algebraic 
soliton at the leading edge. The two bores match at the 
trailing edge of the cnoidal bore, which is defined by ([77| 
with r~ = u l (l — au l ) and r + = u r (l — au r ) 



12u r (l-au r )-6u\l-au l ) 



(98) 



The amplitude at the matching point for the normal 
composite bore is (see ([85]) . (1901) ) 



a* = 2(R 2 (s*) - J?i(s*)) = 4 



1 

~2a~ 



(99) 



The same result obviously holds for the reversed com- 
posite bore, for which we use (|96|) and (|97| . The trailing 
edge speeds for both types of composite bores is given by 



B. Classification for a > 

We can now proceed with the full classification of the 
solutions to the step problem. 
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FIG. 5: (Color online) Parametric map of solutions of the step 
problem for the Gardner equation with a > 0. The resolution 
diagrams corresponding to each of the cases on the plane of 
the initial step parameters u~ and u + are the following: Re- 
gion 1: {u~ UB -> u + }; Region 2: {u~ <- UB («*) SB 
u + }; Region 3: {u~ RW (u*) SB u + }\ Region 4: 
^- RW u + }; Region 5: {u~ <- UB Region 6: 

UB (m*) SB u+}; Region 7: RW (u*) SB 
Region 8: RW — > u + }. In all relevant cases the 
intermediate state u* = 1/q — u + . 



We first present a detailed classification of solutions to 
the evolution of a step problem for the Gardner equa- 
tion with a > 0. The parametric map of solutions is 
constructed on the (u + ,u~) plane of the initial step pa- 
rameters (see Fig. 5). The whole (u~ , w + )-plane is split 
into 8 regions, each corresponding either to one of the 
basic patterns listed in Section IV A ("pure" solutions) 
or to the combination of two patterns ( "composite" so- 
lutions). To represent the result of the evolution of an 
initial step for each region we shall be using symbolical 
diagrams introduced in the previous Section. Say, the 
resolution diagram for Region 7, {u~ UB — >• (u*) <— SB 
u + }, denotes a normal undular bore connecting the left 
state u~ with an intermediate state u* which is further 
connected to the right state u + via the reversed soliborc. 

The lines separating the regions are: 

• u~ = u + — separates the regions of pure undular 
bores and pure rarefaction waves; 

• u~ = l/a— u + — corresponds to the steps resolving 
into single solibores and separates the regions of 
composite solutions of different types: undular bore 
+ solibore and rarefaction wave + soliborc; 

• u~ = 1 /2a — separates regions of pure (undular 
bore) and composite (undular bore + solibore) so- 
lutions. 

We note that the classification for a > is qualita- 



tively similar to that presented in [22( for the Myatta- 
Choi-Camassa (MMC) system describing fully nonlinear 
interfacial dispersive waves in a two-layer fluid. In [22j], 
the analytic method of [2(| was used to obtain the loca- 
tions of undular bore boundaries and the leading solitary 
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wave amplitude. However, the full modulation solutions 
are not available for the MMC system due to complex- 
ity of the corresponding Whitham equations. The the- 
ory presented in this paper has an obvious advantage of 
greater simplicity and universality due to the integrable 
nature of the problem and availability of exact analytic 
solutions. At the same time, in the context of internal 
water waves, the Gardner equation with a > 0, being 
a weakly nonlinear approximation of the MCC system, 
is quantitatively valid only for the waves of sufficiently 
small amplitude. 

The classification is most conveniently performed us- 
ing the function w(u) — u(l — au) defining the mapping 
from the dispersionless Gardner equation (fTTj) to the dis- 
persionless limit of the KdV equation (the Hopf equa- 
tion (|72p (see Section IIIC). We shall illustrate each wave 
pattern in the classification by presenting the analytical 
(modulation theory) solutions along with respective di- 
rect numerical solutions of the Gardner equation. In our 
numerics, equation (fTJ was solved using the method of 
lines (see e.g. [37]]) where the spatial derivatives are dis- 
cretised using second order accurate finite difference ap- 
proximation to reduce the governing partial differential 
equation to a system of ordinary differential equations. 
This system is then solved using the fourth order Rungc- 
Kutta method. 



at the trailing edge and 



Region 1, 



UB 



u+} 



Both values u~ and u + lie in the domain where the 
function w(u) = u(l — au) is monotonically increasing so 
there is one-to one correspondence between the disper- 
sionless limits of the Gardner and the KdV equations. 
This suggests that in the Region 1 initial discontinuity 
can be resolved by a single normal "shallow-water" un- 
dular bore of the KdV type with the bright soliton at the 
leading edge and the linear wavepacket at the trailing 
edge. 

We shall use the travelling wave solution (fT3j) and the 
relationships ([5U)) between Uj's and r^'s to construct the 
desired modulated travelling wave solution for the undu- 
lar bore and show that it indeed provides the required 
matching between u~ and u + . The parameters enter- 
ing the modulation solution ([70)1 . ([75)1 can be expressed 
in terms of the initial step parameters using the rela- 
tionships (|59p. 

It follows for the small-amplitude limit (j2"Tj) of the trav- 
elling wave solution (TT3"|) that the trailing edge of the un- 
dular bore (m ~ 0) propagates against the background 
u = U2 = U3 (this can be inferred directly from the or- 
dinary differential equation (O, where we set 112 — > U3 
(see Fig. 2b). Similarly, for the soliton edge m — > 1 we 
have that the leading bright soliton propagates on the 
background u = u\ = u 2 (see (|1 1T|) and Fig. lb). Thus, 
if the step is resolved by a single undular bore we must 
require 



U2 = U\ = U 



(101) 



at the leading edge. 

Considering the same limits in the relationships 
Eqs. ([55} we have 

1 



«2| m=0 = u 3 \ m=0 = — (1 - y/1 - 4ar 3 ) (102) 



and 



«iL=i = "2| m=1 = 2^(l-Vl-4an) ■ ( 103 ) 

According to (|74[) = r~ and n = r + in the Gurevich- 
Pitaevskii solution. Then, from (|100[) ~ ([103[) we have 



r = u (1 — au ) = w(u ), 
r + = u + (l — au + ) = w(u + ) . 



(104) 



Since for the considered Region 1 we have w(u~) > 
w(u + ), the condition ([75)1 of the Riemann invariant or- 
dering is satisfied, therefore our construction is consistent 
throughout. 

The undular bore occupies the region s~t < x < s + t 
where the edge speeds s ± are obtained from (|77[) and ([75]) 



s = 12w + (l — au + ) — 6m (1 — au ), 
,s + = 2u + (l - au + ) + 4u~(l - an"). 

The width of the undular bore is then 
A = (s+ - s~)t = 10(r~ - r + )t 
= 10(u" - u + )(\ + a(u~ + u + ))t . 



(105) 



(106) 



One can see that the Gardner undular bore is wider than 
its KdV counterpart, for which one has A^v = 10(«_ — 
u + )t for the same initial conditions. As a matter of fact 
the KdV result is reproduced when a = 0. The amplitude 
of the lead soliton in the undular bore is (see (fT8)l . ([59]) 
and QUI)) 



= («3-Wi)L=i = 2(m -■ 



(107) 



This result coincides with the classical KdV formula for 
the lead soliton amplitude. 

The constructed solution is illustrated in Fig. 6 where 
the plot (a) of the analytical (modulation theory) solu- 
tion is presented along with the direct numerical solu- 
tion of the Gardner equation (plot (b)). One can see 
that agreement is very good. The presence in the numer- 
ical plot of an extended small-amplitude oscillatory tail 
stretching behind the trailing edge as defined by the mod- 
ulation theory is a well-known feature of undular bore 
solutions observed in the early comparisons of the KdV 
modulation solutions with numerics (see e.g. [HJ) and 
recently studied in detail in (39| . 



Region 2, u + < ^- < u < 1/a 



{u- 



UB 



U 2 = "3 



(100) (u*)SB^u+}. 
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FIG. 7: (Color online) Finding the intermediate state u = u*. 
(a) Region 2 diagram, it* > u~; (b) Region 3 diagram, it* < 



and 



FIG. 6: (Color online) Evolution of an initial discontinuity for 
the Gardner equation with a = 1. Region 1: {it~UB — m + }). 
The step parameters are u~ = 0.3, it + = 0.1. (a) analytical 
solution in the form of a modulated periodic wave; (b) nu- 



merical solution; the analytically found edges x 



H of 



the undular bore are shown by dashed lines, 
respond to t = 300. 



Both plots cor- 



U3 = M4 = u = u 

a 



Then from (1601) we obtain 



«2| m=0 



"3| m=0 = ^ (1 + VI - 4ar 3 ) 



(109) 



(110) 



and 



Since u + and u~ now lie in different regions of mono- 
tonicity of the function w(u) = it(l — au), the Region 1 
solution in the form of a single normal undular bore is 
not able to provide the necessary continuous matching 
between the given states. Instead, a reversed undular 
bore is generated between u~ and the intermediate state 
u* = (1/a — u + ) > it - . This intermediate state is found 
from the condition w(u*) = w(u + ) and corresponds to 
the required boundary value for the Riemann invariant, 



(1 



t + ) but lies in the same as u + region of 



monotonicity of the mapping function w(u) (see Fig. 7a). 
The further connection between u* and u + is provided 
by a normal solibore. The corresponding wave pattern is 
presented in Fig. 8. 

The relationship between the Riemann invariants {r{\ 
and the undular bore parameters {u{\ is now described 
by formulae ([60|l . The trailing edge (m = 0) of the bore 
propagates on a background u = 113 = 112 = u~ . At 
the leading edge (to = 1) we have a dark soliton ([TO)) 
propagating against the background u = 1*3 = 1*4. Thus, 
for the undular bore we have 



"3 1 



U 4 



=0 = ^ (1 + VT= 4anJ 



(111) 



Again, since n = r+ and r% = r_ in the modulation 
solution for the undular bore, we have from (|108[) - (|llip : 



r3 = r = u (1 — au ) = w(u ), 
ri = r + = u*(l — au*) = u + (l — au + ) . 



(112) 



Thus, the expressions for the undular bore speeds remain 
the same (cf. (JTUSJ)): 



s = 12k (1 — au ) — 6u (1 — au ) . 
s + = 2u + (l - au + ) + 4u~(l - air) . 



(113) 



The front solibore connecting the states u* and u + is 
described by formula ([2TJ]). The solibore speed is given 
by 



s k = — + 2au + u* = — + 2u + (l - au A 
a a 

Now it is not difficult to sec that 



(114) 



U3 = u 2 



(108) 



s k - s+ 



Au (u 



2a 



) 2 >0, 



(115) 
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FIG. 8: (Color online) Evolution of an initial discontinuity 
for the Gardner equation with a = 1. Region 2, {u~ <— 
UB («*) SB —J- u + }. The initial step parameters are 
u~ — 0.7, u + = 0.1. (a) analytical solution in the form 
of a modulated periodic wave connected to a solibore; (b) 
numerical solution of the Gardner equation; the analytically 
found edges = s t and x k = s k t of the undular bore and 



solibore respectively are shown by dashed lines, 
correspond to t = 300. 



Both plots 



that is the solibore always propagates ahead of the undu- 
lar bore as expected. Since there is a qualitative change 

of the wave pattern in the transition from Region 1 to 
Region 2 it is necessary to look closer at what happens 
near the boundary between these two regions determined 
by the value of the left state u~ = l/2a. The change of 
the pattern is illustrated in Fig. 9 where several numer- 
ical solutions of the Gardner equation with a — 1 arc 
presented for the evolution of initial discontinuities with 
the same right state, u + = 0.1, while the left state it- 
was taken in the range u_ = 0.49 < 1/2q (Region 1) to 
U- = 0.58 > l/2a (Region 2). 

One can see from (|115[) that for u~ = l/2a the speed of 
the solibore coincides with the speed of the leading soli- 
ton in the undular bore so the "borderline" wave pattern 
in Fig. 9b can be interpreted in both ways: as a normal, 
"bright" , undular bore or, equivalently, as a reversed, 
"dark", undular bore with an attached solibore. When 
we increase u~ , the solibore separates from the undular 
bore, which in its turn acquires the reversed waveform 
with a distinct dark soliton structure near the leading 
edge (see Fig. 9d). 




100 150 200 250 300 350 400 450 



FIG. 9: (Color online) The transformation of the undular bore 
structure from the normal, "bright soliton" , pattern in Region 

1 (plot a) to the reversed, "dark soliton", pattern in Region 

2 (plot d). Numerical simulations of the Gardner equation 
with a — 1. The downstream state, u + = 0.1, is the same 
for all cases, the upstream state u~ is taken in the range 
u~ = 0.49 < l/2a (Region 1) to u~ = 0.58 > l/2a (Region 
2), t=300 



Region 3, 1/a—u < u + < ^ < u , u +u + > 1/a. 
{u~ <- RW (it*) SB -> u+} 

This region is analogous to Region 2, since the val- 
ues u~ and u + again lie in different domains of mono- 
tonicity of the function w(u), thus a single- wave reso- 
lution is not possible. However, now the intermediate 
state u* = 1/a — u + < u + (see Fig. 7b) so a reversed 
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FIG. 10: (Color online) Evolution of an initial discontinuity 
for the Gardner equation with a = 1. Region 3: {u~ 
RW (it*) SB — !> u + }. The initial step parameters are: vT = 
1.0, u + = 0.3; (a) analytical (dispersionless limit) solution; 
(b) numerical solution of the Gardner equation. Dashed lines 
in (b) correspond to the analytically found locations of the 
rarefaction wave edges x l ' r = s l ' r t and solibore position x k = 
s t. Both plots are made for t = 1500. 
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FIG. 11: (Color online) Evolution of an initial discontinuity 
for the Gardner equation with a = 1. Region 4: {u~ RW 
u + }. The initial step parameters are: u~ = 0.9, u + = 0.6. 
(a) analytical (dispersionless limit) solution; (b) numerical 
solution. Dashed lines in (b) correspond to the analytically 
found locations x Lr of the rarefaction wave edges. Both plots 
correspond to t = 600. 



rarefaction wave is generated instead of reversed undu- 
lar bore. The solution for the rarefaction wave is given 
by formula (|81|) . where ui = u~ and u r ~ u* . The 
solibore solution connecting u* and u + is the same as 
in Region 2. The analytical and numerical plots cor- 
responding to Region 3 are presented in Fig. 10. We 
note that, since maxit + (l — au + ) = j- , the speed of 
the solibore Sk = — + 2u + (l — au + ) is always greater 
than that of the right edge of the rarefaction wave, 
s + = s r = 6u+(l - au + ) (see ([82])). At the bound- 
ary between Regions 3 and 4, when u + = l/2a, we have 
Sk = Sr and the solibore gets "attached" to the right edge 
of the rarefaction wave. 

Region 4, 1/2q < u + < u~ . {ir ^- RW u + }. 

A single reversed rarefaction wave is produced. It is 
described by the solution (|81j) with ui = u~ , u r = u + . 
The corresponding typical analytical and numerical solu- 
tion are presented in Fig. 11. Note the small-amplitude 
dispersive wave train seen to the left of the rarefaction 
wave in the numerical plot in Fig. lib. This wavetrain is 
necessary to resolve weak discontinuity at the left edge of 



the rarefaction wave, while another weak discontinuity at 
the right edge is smoothed out (see [l?} for the detailed 
description of similar effects in the KdV theory). 

Region 5-8. The wave patterns corresponding to 
Regions 5-8 represent the "reflections" of the patterns 
arising in Regions 1-4. More precisely, the counter- 
part solutions correspond to the "opposite" regions in 
the parametric map in Fig. 5 and are related to each 
other by the transformation ([5]). From this viewpoint, 
the solutions for Regions 5-8 are not "new". At the 
same time, the change of the solution polarity (e.g. from 
the "bright undular bore" to the "dark undular bore") 
due to the change of initial data is not trivial physically 
so it deserves separate description. For this reason and 
for the reader to be able to identify the arising wave pat- 
terns directly, without the need to invoke intermediate 
transformations, we shall proceed with the descriptions 
of the Regions 5 - 8 in the same format that was used for 
Regions 1-4. 

Region 5, l/2a < u~ < u + < i. {u~ <- UB u + } 
In this region both values u~ and u + are in the domain 
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FIG. 12: (Color online) Evolution of an initial discontinuity 
for the Gardner equation with a = 1. Region 5: {u~ <— UB 
u + }. The initial step parameters are: u~ — 0.6, u + — 0.9. (a) 
analytical (modulation theory) solution; (b) numerical solu- 
tion of the Gardner equation. Dashed lines in (b) correspond 
to the analytically found locations x of the undular bore 
edges. Both plots correspond to t — 300. 



FIG. 13: (Color online) Evolution of an initial discontinuity 
for the Gardner equation with a = 1. The initial step param- 
eters are: u~ = 0.3, u + = 0.9. Region 6: {u~ UB — > (w*) <— 
SB u + }). (a) analytical (modulation theory) solution; (b) nu- 
merical solution of the Gardner equation. Dashed lines in (b) 
correspond to the analytically found locations of the undular 
bore edges x and the solibore x k . Both plots correspond to 
t = 200. 



when the function w(u) = u(l — cm) is monotonically de- 
creasing so a single reversed undular bore is produced. 
The modulation description of this undular bore is iden- 
tical to that in Region 2 but to obtain the oscillatory 
structure one now needs to use the relations (|60|) between 
the Riemann invariants and parameters of the periodic 
solution (|13p . The amplitude of the lead dark soliton 
a + = («4 — U2)|m=i is given by by the same expression 
(|107p . The analytical and numerical solution plots for 
Region 5 are presented in Fig. 12. 

Region 6, ^-u + < u~ < i < u + . {u~ UB -> (u*) 
<- SB u+} 

A combination of the normal undular bore and a re- 
versed solibore is produced. 

The undular bore connects the state u~ and an inter- 
mediate state u* = 1/a — u + and is described by the 
same set of formulae as a single normal undular bore in 
Region 1. The reversed solibore further connects the in- 
termediate state u* with the downstream state u + . It 
is described by formula (|2T)|) with "+" sign. The plots 
of the analytical and numerical solutions are shown in 
Fig. 13. 



Region 7, u <j^<u + <l/a-u. {u RW 
-> (u*) <- SB u+} 

The resolution pattern is similar to that in Region 6 
but, instead of the normal undular bore, a normal rar- 
efaction wave described by formula (|BT))) is generated. 
The boundary states are ui = u~ , u r — u* — 1/a — u + . 
The corresponding analytical and numerical plots are 
presented in Fig. 14. Note that, similar to Region 3, 
the solibore always propagates ahead of the rarefaction 
wave and gets attached to the right edge of the rarefac- 
tion wave when u + = I /2a. 



Region 8, u < u + < 

A single normal rarefaction wave is produced, {u~ RW 
— > u + }. It is described by the solution (|80]l with u l = u~ , 
u r = u + . The plots of the analytical and and numerical 
solutions for Region 8 are shown in Fig. 15. 
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FIG. 14: (Color online) Evolution of an initial discontinuity 
for the Gardner equation with a, = 1. Region 7: {it - RW 
— > (u*) SB u + }. The initial step parameters are: \T — 
0.1, u + = 0.7. (a) Analytical (dispersionless limit) solution 
in the form of a rarefaction wave connected to a solibore; (b) 
Numerical solution of the Gardner equation. Dashed lines 
in (b) correspond to the analytically found positions of the 
solibore x k and rarefaction wave boundaries x l ' r Both plots 
correspond to t = 700. 



C. Classification for a < 

Now we present the parametric map of solutions of 
the step problem for the Gardner equation with a < 0. 
The most significant change in the structure of solutions 
compared to the case a > is that the composite so- 
lutions now contain trigonometric undular bores rather 
than solibores. The plane (u + , u~) of the initial step pa- 
rameters is again split into 8 regions (see Fig. [To]) . The 
lines separating different regions are: 

• u~ = u + separates the regions of pure undular 
bores and pure rarefaction waves; 

• u~ = 1/a — u + corresponds to the steps resolving 
into single trigonometric bore solutions and sepa- 
rates the regions of composite solutions of different 
types: undular bore + trigonometric undular bore 
and rarefaction wave + trigonometric bore; 

• u~ = l/2a separates regions of pure (undular bore) 
and composite (undular bore + trigonometric bore) 



FIG. 15: (Color online) Evolution of an initial discontinuity 
for the Gardner equation with a = 1. Region 8: {u~ RW — s- 
u + }. The initial step parameters are: u~ = 0.1, u + = 0.4. 
(a) Analytical (dispersionless limit) solution in the form of 
a rarefaction wave; (b) Numerical solution of the Gardner 
equation. Dashed lines in (b) correspond to the analytically 
found positions x l ' r of the rarefaction wave boundaries. Both 
plots correspond to t = 500. 



solutions. 

Let us now describe in some detail the wave structures 
corresponding to different regions in Fig. 1161 

Region 1. u + < u~ < {u~ «- RW u + }. 

Both values u~ and u + lie in the domain where the 
function w(u) = u(l — au) is monotonically decreasing 
so there is one-to one correspondence between the disper- 
sionless limits of the Gardner and the KdV equations. 
This suggests that the Region 1 initial discontinuity is 
resolved by a single "reversed" simple rarefaction wave. 
The rarefaction wave is described by the solution (|5Tj) 
with u l = vT and u T = u + . A typical solution for Re- 
gion 1 is shown in Fig. ??. The plot for the corresponding 
analytical solutions is almost identical to that shown in 
Fig. ?? with the exception for the small-amplitude wave 
train at the left corner of the rarefaction wave so we do 
not present it here. The analytically found boundaries of 
the rarefaction wave are shown by the dashed lines. 



Region 2. u + 

(u*) <- RW u+} 



< 



l 

2n 



< U~ < 



{u~ TB 



20 



\j 3 


u , 

4 




2 




5 

1 


1 / 

/ 8 

/ 2" 


7 

If 


2a 

6 

\ 



FIG. 16: (Color online) Parametric map of solutions of the 
step problem for the Gardner equation with a < 0. The res- 
olution diagrams corresponding to each of the cases on the 
plane of the initial step parameters u~ and u + are the fol- 
lowing: Region 1: {u~ <— RW tt + }; Region 2: {u~ TB — > 
(u*) <- RW u + }; Region 3: {u^ (TB|UB) -> Re- 
gion 4: {tt~ UB — > it + }; Region 5: RW — > 
Region 6: {«" <- TB («*) RW -> n+}; Region 7: 
{iT <- (TB]UB) Region 8: {u~ <- UB «+}. In 
all relevant cases the intermediate state u* = 1/a — u~ . 




(a) 

-loo -100 100 200 300 400 500 
X 

















r 




(b) x 1 


X 





f 00 -100 100 200 300 400 500 
X 



FIG. 17: (Color online) Evolution of an initial discontinuity 
in the Gardner equation with a — — 1: Region 1: {u~ -s— RW 
u + }. The initial step parameters are: u~ = — 1, u + = —2.5. 
(a) Analytical (dispersionless limit) solution; (b) Numerical 
solution of the Gardner equation. Dashed lines in (b) cor- 
respond to the locations of the rarefaction wave edges found 
analytically. Both plots are made for t = 20. 
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FIG. 18: (Color online) Evolution of an initial discontinuity 
for the Gardner equation with a = — 1 for t = 3. Region 2. 
{u~ TB — > (u*) <— RW u + }. The initial step parameters: 
■uT = 1, u + — —3. (a) Riemann invariants Ri, R2 and R3. 
(b) Analytical (modulation theory) solution, (c) Numerical 
solution of the Gardner equation. Dashed lines on the nu- 
merical plot correspond to the analytically found boundaries 
between different parts of the wave pattern. 



Since u + and u~ now lie in different regions of mono- 
tonicity of the function w(u) = u(l — cm), the Region 1 
solution in the form of a single reversed rarefaction wave 
is not able to provide the necessary continuous matching 
between the given states. Instead, one needs to introduce 
a trigonometric undular bore joining the left constant 
state u~ with rarefaction wave at the level u* = l/a — u~ 
(see Fig. rrS)) . The intermediate state u* is found from the 
condition w(u*) = w(u~) (cf. condition w(u*) = w(u + ) 
for a > 0). The modulation description of the relevant 
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FIG. 19: (Color online) Evolution of an initial discontinu- 
ity for the Gardner equation with a — —2. Region 3: 
{u~ (TB|UB) — > u + }. The initial step parameters are: 
u~ = 0.2, u + — —0.5. The plots are shown for t — 50. (a): 
Riemann invariants Ri, Ri and R3. (b): Analytical (modula- 
tion theory) solution, (c): Numerical solution of the Gardner 
equation. Dashed lines on the numerical plot correspond to 
the analytically found boundaries between different parts of 
the wave pattern. 



(normal) trigonometric bore was constructed in Section 
IVA in terms of the Riemann invariants R3 > R2 > Ri 



The speeds s and s* of the trailing and the leading 
edges of the trigonometric bore, and the speed of the 



leading edge of the rarefaction wave s + are: 
3 



— 6it (1 — au ), 



s* = 6u*(l — au*) = 6u (1 — au ), 
s + = 6u + (l — au + ). 



(116) 



The analytical and numerical solutions along with the 
plot for the Riemann invariants Ri, R2, i?3, are shown in 

Fig.ua 

At u~ = 1/a — u + > u + the rarefaction wave disap- 
pears and one obtains a single normal trigonometric bore 
as a result of the step evolution The relevant analytical 
description was presented in Section IVA. 



Region 3. —— u 



< u 



< {„- (TB|UB) 



->• u 



In the Region 3 we get a composite undular bore con- 
sisting of normal trigonometric and cnoidal parts match- 
ing at the point of the trailing edge x* = s*t of the 
cnoidal bore. The modulation solution for the entire com- 
posite bore is conveniently described in terms of the Rie- 
mann invariants R3 > R2 > Ri (|83|) . The corresponding 
analytical and numerical plots are shown in Fig. [TH1 

The characteristic speeds for this region are 



3 

a 



— 6u (1 — au ), 

s* = 12u+(l - au + ) - 6u~(l - au~), 
,s + = 2u + (l - au + ) +4u - (l - au~). 



(117) 



We note that the boundary x* = s*t between the 
trigonometric bore and the cnoidal bore parts in the com- 
posite bore solution can be naturally defined only in the 
framework of the averaged (Whitham) equations, where 
it represents a characteristic separating two regions with 
qualitatively different behaviour of the modulation solu- 
tion (the modulation solution has a weak discontinuity at 
x = x* - see Fig. [T9b). Due to the asymptotic nature of 
the modulation equations (the phase is washed out), this 
separating line cannot be consistently identified on the 
level of the genuine (rapidly oscillating) solution of the 
governing equation. The situation here is similar to that 
with the definition of the trailing edge of a standard un- 
dular bore: the trailing edge cannot be identified with a 
particular point in the bore but is rather associated with 
the linear group velocity characteristic of the modulation 
equations, hence the already mentioned noticeable dif- 
ference in the behaviour of the asymptotic (modulation 
theory) solution and that of the full numerical solution 
near the trailing edge, (cf Figs. [T5b and ITT)b) . 

Region 4. u~ > u + > {u~ UB ->• w+}. 

Both values u" > u + lie in the domain where the func- 
tion w(u) = u(l — au) increases, so the resolution occurs 
via a single normal cnoidal undular bore; see Fig. 1201 To 
be consistent with other plots in this section we present 
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FIG. 20: (Color online) Evolution of an initial discontinuity 
for the Gardner equation with a = — 1. Region 4: {u~ UB 
— > u + }. The initial step parameters are u~ = 1.5, it + = 0.5. 
The plots are shown for t = 3. (a) Riemann invariants Ri, 
R2 and R3. (b) Analytical (modulation theory) solution, (c) 
Numerical solution of the Gardner equation. Dashed lines 
on the numerical plot correspond to the analytically found 
boundaries between different parts of the wave pattern. 



the Gurevich-Pitaevskii solution (j74j) . (|75|) for the undu- 
lar bore in terms of the Riemann invariants {Rj} rather 
than {rj}. The one-to-one correspondence between these 
two sets of the Riemann invariants in Region 4 is given 
by relations (|83|) . 

The edge speeds are (see f77]). ([75]) ): 



s = 12u (1 — au ) — 6it (1 — au 
s + = 2u~(l - au~) + 4u + (l - au + ). 



(118) 



Regions 5-8. Similar to the classification for a > 
described in the previous section, the solutions for Re- 
gions 5-8 can be obtained by applying the transforma- 
tion ([5]) to their counterparts from the "opposite" regions 
of the parametric map in Fig. 16. Again, for the con- 
venience of identification we present them below in the 
same format. 




400 600 
X 



2.5 



0.5 



(b) 






1 


r 




X 


X 





200 



400 600 
X 



800 



1000 



FIG. 21: (Color online) Evolution of an initial discontinuity 
in the Gardner equation with a — —1. Region 5: {u~ RW 
— > u + }. The initial discontinuity parameters are: u~ = 0.5, 
u + = 2; (a) Analytical (dispersionless limit) solution; (b) Nu- 
merical solution of the Gardner equation. Dashed lines cor- 
respond to the analytically found locations of the rarefaction 
wave edges. The plots are shown for t = 20. 



Region 5. 



77- < u < u + . 

2a 



{w RW -> u+}. 



Both values u~ < u + lie now in the domain where the 
function w(u) = u(l — au) decreases, so the resolution 
occurs via a single normal rarefaction wave described by 
([50]) with u l =u~, u r = u+. 

The edge speeds are 



6it (1 — au ), 



6u+(l-cra + ). (119) 



The corresponding numerical solution is shown in Fig. [5T] 
along with the boundaries of the analytical RW solution 
marked by dashed lines. 



Region 6. — — u + < 



TB (u*) RW 



"}• 



< 2^ 



{> 



The resolution pattern corresponding to this region is 
similar to that in Region 2 but now the resolution occurs 
via the combination of the reversed trigonometric bore 
and normal rarefaction wave joined at the level at the 
level u* = l/a — u~ . The modulation description of 
such a bore is constructed in Section IVA in terms of the 
Riemann invariants R\ > R2 > R3 (|84p. The relevant 
plots are shown in Fig. [5U 
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FIG. 22: (Color online) Evolution of an initial discontinu- 
ity for the Gardner equation with a = — 1. The initial 
step parameters: u~ = —1, u + = 0.4 (Region 6. {it - <— 
TB («*) RW — > u + })- (a) Riemann invariants Ri, R2 and 
i?3 . (b) Analytical (modulation theory) solution, (c) Numer- 
ical solution of the Gardner equation, t — 50. Dashed lines 
on the numerical plot correspond to the analytically found 
boundaries between different parts of the wave pattern. 



FIG. 23: (Color online) Evolution of an initial discontinu- 
ity for the Gardner equation with a = —2. Region 7. 
{u~ (TB|UB) u + }). The initial step parameters are: 
u~ = —1.6, u + — 0.6. The plots are shown for t — 2. (a): 
Riemann invariants Ri, R2 and R3. (b): Analytical (modula- 
tion theory) solution; (c): Numerical solution of the Gardner 
equation. Dashed lines on the numerical plot correspond to 
the analytically found boundaries between different parts of 
the wave pattern. 



The edge speeds are equal to 



s = 6u (1 — au ), 

a 

s* = 6w~(l — au~), 
s + = 6u + (l — au + ). 



(120) 



The part of the line u~ — 1/a — u + , where u~ < u + , 
separating regions 6 and 7 corresponds to a pure reversed 
trigonometric bore described in Section IVA. 



Region 7. 

(TB|UB) u+} 



< 



1 

2d 



< M" 1 



< 



Region 7 corresponds to the formation of a composite 
reversed trigonometric-undular bore (cf. Region 3 for the 
counterpart normal resolution pattern) . The correspond- 
ing plots are shown in Fig. 1231 
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This region corresponds to the to formation of a re- 
versed cnoidal bore (cf. Region 4); the plots are pre- 
sented in Fig. [24j Note that the Riemann invariants 
Ri > R2 > R3 were used in the construction of the 
modulation solution for the reversed bore; as a matter 
of fact, it is equivalent to the Gurevich-Pitacvskii solu- 
tion ([71)) . (|75|) in the original variables > r% > r\ (sec 
Fig. 3). 

The edge speeds are 
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FIG. 24: (Color online) Evolution of an initial discontinuity 
for the Gardner equation with a = —1. Region 8. {u~ 
UB u + }. The initial step parameters are u~ = —1.5, u + — 
—0.7. The plots are shown for t = 10. (a)Riemann invariants 
Ri, R2 and R3. (b) Analytical (modulation theory) solution, 
(c) Numerical solution of the Gardner equation. Dashed lines 
on the numerical plot correspond to the analytically found 
boundaries between different parts of the wave pattern. 



The edge speeds are 



s = 6u (1 — au 

a 

c* — 10,,+ 



12u + (l — au ) — 6u (1 — au ), 
2u + (l - au + ) + 4u~(l - em - )- 



(121) 



Region 8. u < - - u + , u + < {u <-UBu + l. 



12u + (l — au + ) — Qu (1 — au ), 

x x _ _ (122) 

2it + (l - + Au (I - au ). 



V. CONCLUSIONS AND OUTLOOK 

We have constructed a full analytical description of 
the step problem for the Gardner equation ((T|) for both 
signs of the coefficient a before the cubic term. The 
complete classification of arising solutions for different 
parameters u + ,u~ defining the initial step ([B]) includes 
16 possible cases (8 for each sign of a) . Each sector on the 
(u + , u~) plane of the parametric map of solutions corre- 
sponds to a unique wave pattern representing one of the 
following: undular bore, a rarefaction wave, a solibore, a 
trigonometric bore; or a combination of two of the above 
wave structures. The wave pattern arising in each case 
depends on the position of the initial step parameters 
u~,u + relative to each other and to the turning point 
l/(2a) of the function 6u(l — au) defining the character- 
istic speed of the dispersionless limit of the Gardner equa- 
tion. The analytical description of undular bores is made 
using the Whitham modulation theory. The observed 
rich phenomenology of solutions arising in the step prob- 
lem for the Gardner equation is due to the fact that the 
modulation Whitham system associated with the Gard- 
ner equation, unlike that for the KdV equation, is neither 
strictly hyperbolic nor genuinely nonlinear. Our analyt- 
ical solutions are supported by numerical simulations. 

One of the important applications of the obtained so- 
lutions is an analytical description of transcritical flow 
in a stratified fluid in the framework of the forced Gard- 
ner equation (cf. [2^, Hoj]). Other possible applications 
include the consideration of the interaction of internal 
undular bores with variable topographies (cf. [4(|) and 
the description of the perturbed modulation regimes for 
internal waves (e.g. due to the inclusion of weak dissi- 
pation). In the latter case, the description will require 
a perturbed modulation theory approach developed in 
4l| . The obtained classification will also provide a guid- 
ance for the similar classifications for fully nonlinear non- 
integrable counterparts of the Gardner equation (such as 
Myatta-Choi-Camassa system (3f|) to which the analytic 
technique of the undular bore description developed in 
[20j is applicable (see (22j). 
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